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Abstract. The subject of the paper is the role of the 
massive stellar cluster in the activity phenomenon and in 
the structure of active galactic nuclei. We introduce a sim- 
ple model of stellar dynamics in the internal part of the 
cluster, which allows us to include both the star-disk and 
the star-star interactions. It is shown that the properties 
of the distribution of stars in the vicinity of the black 
hole are determined both by the interaction of the stars 
with the accretion disk and by the pair gravitational and 
contact interaction between the stars. We calculate the 
distribution of stars in the central parts of the cluster and 
we discuss possible effects of stellar mass-loss due to the 
star-disk interaction. Finally, we study the implications of 
the central cluster for active galactic nuclei activity. We 
model the broad line region assuming that the gaseous 
wakes, following stars after each disk crossing, play the 
role of the broad line region clouds, and we calculate the 
corresponding line profiles. We also analyze the contribu- 
tion of star-star and star-disk collisions to active galactic 
nuclei variability. 

Key words: Galaxies: active. Galaxies: Seyfert, Stellar 
dynamics. Ultraviolet: galaxies. X-rays: galaxies 



1. Introduction 

The idea of an important role of the compact stellar clus- 
ter (hereafter CSC) in creation and evolution of active 
galactic nuclei (hereafter AGN) was formulated at the very 
early stage of AGN investigation (e.g., Shklovskiy 1964), 
and the evolution of a gas sphere containing a CSC was 
first considered in the work of Unno (1971). The dissi- 
pative interaction of the compact stellar cluster with the 
massive central object in the form of a super-star was 
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considered by Vilkoviskij (1975) and Hara (1978), and the 
tidal interaction of stars with the central object in the 
form of a massive black hole (hereafter MBH) - by Hills 
(1975). The evolution of the dense, non-rotating stellar 
cluster was discussed by Saslaw (1966), Bisnovatyi-Kogan 
(1978), and others. The interaction of stars with the ac- 
cretion disk (hereafter AD) was considered in the works 
by Vilkoviskij and Bekbosarov (1981), Vilkoviskij (1983), 
Zentsova (1983), Syer et al. (1991). More detailed investi- 
gations of the evolution of the stellar orbits, crossing the 
AD were presented in Artymowicz et al., (1993), Karas 
and Vokrouhlicky (1993), Ranch (1995), Vokrouhhcky & 
Karas (1993, 1998), Subr & Karas (1999), and others. The 
role of the outer parts of AD (which can be gravitation- 
ally unstable) in the global evolution of AGNs through 
efficient star formation was discussed by Collin & Zahn 
(1999). 

A general investigation of the problems of CSC struc- 
ture and evolution, taking into account the interaction of 
the star cluster with the accretion disk, demands calcula- 
tions of very complicated numerical models. In contrast, 
the main aim of the present paper is a preliminary anal- 
ysis at a qualitative level, but including both star-star 
and star-disk interactions. The star-disk and the star-star 
interactions were not considered before as acting simulta- 
neously. Ranch (1995, 1999) considered and numerically 
calculated the cases when either star-disk or star-star col- 
lisions dominate; we consider here a model which is very 
simple but includes the joint action of the star-disk and 
star-star interactions. This permits us to find the space 
distribution of the stars around ADs. As an application of 
this result, we consider star-related model of Broad Emis- 
sion Line Regions (hereafter BELRs) and the possible role 
of the inelastic interactions of stars both with the disk and 
other stars in the variability of AGNs. 

The plan of the paper is as follows. In Sect. 2 we con- 
sider the stellar dynamics in the central part of CSC in 
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AGNs. In Sect. 3 wc derive the strueture of the stellar 
cluster in the AD region, including the role of contact stel- 
lar collisions and some effects of the mass loss of the stars. 
The application of the results to modeling BELR and vari- 
ability are discussed in Sects. 4 and 5, correspondingly. 
The discussion of the results and conclusions are given in 
Sect. 6. 

2. Stellar dynamics in the central part of CSC 

2.1. The star- disk interaction 

We assume that there are massive and compact stellar 
clusters in the centers of AGNs. Observations of the near- 
est galaxies show clearly that dense stellar clusters do exist 
in their centers, with sizes close to a parsec, with masses 
10^ - 10* Mq, and central stellar densities 10^ - 10^ pc'^. 
It is reasonable to suppose that in AGNs the density in 
CSCs is increased, and in the case of the more power- 
ful AGNs the central stellar density can reach values of 
up to 10* — 10^ pc~^, as follows in particular from the 
gas-dynamic solutions for the broad absorption line QSOs 
(hereafter BAL QSO, see Vilkoviskij et al. 1999). In con- 
trast to the case of close galaxies with low activity, the ob- 
servational data about CSCs in luminous AGNs are very 
poor; we will assume that those CSCs are massive, rotate, 
and consist predominantly of low-mass stars. 

For the sake of qualitative analysis, wc introduce sev- 
eral simplifications in our consideration: (1) we consider 
equal stellar masses, (2) we assume that all stellar orbits 
are almost circular, which is a good approximation in the 
internal part of the CSC only (see below), (3) we intro- 
duce some effective values of the main physical parameters 
- the inclination angle of stellar orbits to the plane of AD, 
and the stellar density (see the next subsection), and we 
investigate the radial dependence of these parameters. 

The star-disk (hereafter also s-d) interactions were 
considered in a series of studies. The main, and quite obvi- 
ous result of s-d interaction is that the orbits of stars that 
have been captured into bound orbits crossing the accre- 
tion disk evolve from eccentric to almost circular, with 
diminishing radius and inclination angle (e.g., Syer et al. 
1991). The planes of the orbits approach the plane of the 
AD in the process of the orbit evolution. As a result, the 
stars on these evolved orbits move mostly in the direc- 
tion of the AD rotation, due to the dissipative interaction 
of the stars with the disk. For a single star, this process 
leads finally to the incorporation of the stellar orbit into 
the disk plane, and the further fate of such stars is un- 
defined (Artymowicz et al. 1993). In the next sections we 
will show that consideration of the collective star inflow 
through the disk instead of the single stellar orbit evo- 
lution change the result: due to high stellar density, the 
star-star interaction scatters stellar orbits, thus counter- 
acting the capture of stars by the disk. This effect allows 
the CSC to reach some equilibrium, with fluctuating but 



non- vanishing inclination angle of the orbits, and the dis- 
sipative inflow of the stars continues to the center of the 
AD. 

In the present section we will shortly summarize the 
properties of the individual star-disk interaction. For sim- 
plicity we will take into account two above-mentioned ba- 
sic theoretically deduced properties of the orbits, crossing 
the disk: we assume that all orbits are close to the circular 
shape, and that the inclination angles of the orbits to the 
disk are small. 

Let us consider a star in an orbit, crossing the AD in 
the zone where the gravity of the MBH dominates. The 
velocity of a star, crossing the AD at a distance R is 

V; = qc{RJRY'^ = gc(10-3m8/r_2)^/^ (1) 

where c is the velocity of light, i?g = 2GMbh/c^ is the 
gravitation radius of MBH, mg = Mbh/IO^M© is the 
MBH mass, r_2 = i?/(0.01 pc) is the distance from the 
MBH in units of 0.01 pc, and 1 > g > l/y/2 is the co- 
efficient for the Keplerian orbits: q is close to 1 for the 
very stretched orbits, and q = 1/^/2 for the circular ones. 
In our approximation of almost circular orbits we take 
q = l/-\/2, and ~ ^JGM/R. So in convenient units 
v,^ = y,/(109 cm s-i), 

v^c, ~ 0.67(m8A-2)'/'. (2) 

The normal and tangential velocity components with 
respect to the thin disk are = 14 sini and Vt = Vs cost, 
where i is the orbit-to-disk inclination angle. The gas mo- 
tion in the disk is also close to the circular one, Vd — 14; 
so, the relative velocity of the star and gas in the Kep- 
lerian disk is K-o — + (Vd — 14)^]^^^, and wc obtain 
Wre9 — Vs%2^/'^{1 — cosi)^/^. Note that this velocity is su- 
personic in the disk with ~ 10^ K even if the inclination 
angle is very small, down to i ^ 0.1°. 

It is obvious that the drag force acting at the star in 
the disk is proportional to the cross-section of the star, 
the disk density, and the square of the relative star-disk 
velocity, but the coefficient of the proportionality is known 
with an accuracy of a factor of a few; hereafter we will de- 
note it with Q. Several attempts to estimate Q (Zentsova 
1983, Ostriker 1983, Zurek et al.l994) gave slightly dif- 
ferent results. We will use for the energy loss of the star 
crossing thin AD the expression 

A£ = Q7ri?2Sd(ys)V(^), (3) 

where R^ is the radius of the star, is the surface 

mass density of the AD, and (t>{i) is the function of the 
inclination angle. Within our approximations, one has 
^ii) ~ [{Vn? - {Vd - Vtf]/{V,f = 2cosi(l - cosi) = 
2sin?^ cosf/(l -I- cosi). 

Physically this means that the star gives its velocity 
(in the system co-rotating with disk) to the column of the 
disk matter with cross-section close to that of the star, 
and with height equal to the disk thickness. This brakes 
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the normal component of the star velocity, but accelerates 
the star in the direction of the disk rotation. 
The energy generated by the drag force is 

£^diss ~ QiTRl^i{V,,f = Q7ri?2Sd(V;)22(l - cosz) (4) 

which differs from /S.E with 1/ cosi; we will denote 2(1 — 
cosi) = (t>i{i) = cosi. 

Hereafter we will use the Solar units Ts = Rs/Rq', 
nia = Ms/Mq; so, the cross section of a star is wR^ ~ 
1.54 • lO^^rg cm^. We also define the usual values mg = 
Mbh/1O^M0, m = {dM/dt)/MEd, Mcd ^ rngMeyr-i. 

Supposing that the stellar orbits are close to circu- 
lar, we can approximate the time interval between two 
crossings of the disk by a star as AT = nR/Vg 1.4 • 

inS -1/2 3/2 

lO^mg r_2 s. 

Denoting = 10''S4 g cm~^, one can define the en- 
ergy dissipated per crossing 

AE ~ 0.7 • lO^^Qr^.^ifnsrZlMi) erg, (5) 

and the averaged rate of the energy change of a star, 
AE/AT is 

dE/dt ~ 0.5 • lO^^Qr'^Eiml^'^rZl^^Hi) erg s'^. (6) 

The commonly accepted models for the AD are so- 
called a-disks, the Shakura-Sunyaev model (Shakura & 
Sunyaev 1973, hereafter SS), and the Novikov & Thorne 
(1973) model (hereafter NT); the last was used in Rauch 
(1995). On the other hand, physical assumptions of the 
models, the geometrically thin and optically thick of the 
a-disks theory, are not universal. Other processes like en- 
ergy advection, temperature instability, self-gravity and 
the CSC gravity are essential in different cases (see the 
recent review of Park & Ostriker, 2001). 

We will use mostly the NT model, as it contains a 
simple radial dependence of the surface density of the AD; 
according to Rauch (1995) the surface density of the AD 
is 

S4 ^ 2.7 • 102m^/'a:^/'(m_l)VlOr:^/^ (7) 
Then the energy dissipation per crossing is 

~ 1.9 • 10^'^Qr'',ml^^aZTm_/l\zl^Ui{^) erg, (8) 
and the average rate of the energy change of a star is 

dE/dt ~ 1.35 • 10^^Qr'^my^rZl^^^(l>(i) erg s'K (9) 

According to the SS model, the parameters of the 
accretion disk behave differently in three regions (a, b, 
c), and we have from SS the positions of a/b and b/c 
borders: rab-2 = 2.6mf /''aV2i(^^^^,)i6/2i. = 

18.9m8(rni^bc)^^^; where rab-2 = -Rab/lO"^ pc, cp = 
l-r-i/2. 



The surface density and the height scale in the three 
zones may be written as 

a)5] ~ 2.8 • 10''(a m ipy^m^^^^r^J^ S cm-^ 

~ 3.2 • lO" • mms^ cm, (10) 

6)1] ~ 2.1 -1050-4/5^4/5^^ ^^3/5^-3/5^ 

z c 0.85.10i4a-i/iOmf/2°(m)2/5^iAV2V20^ 

c)S ~ 3.11 • 105a-4/5(m ^)y'°mf^\f/\ (11) 
z ~ 0.67-10"a-i/5m-^/^°(TO(p)3/20r-!/| 

With these relations, the energy dissipation per cross- 
ing in the a, b,and c zones of the disk are (in lO''^ erg): 

a)AE ~ 0.2Qr'^{a my^ms^^^rl/^M^), 
^))A£~ 1.47Qr2a-^/5(m)3/5m^/V:f^</.i(i) (12) 
c)AE ~ 2.mQr^,a-^/\my/^^mf/^\zl^Ui{i)- 

The average energy losses of a star for the three zones 
are (in lO^^erg s~^): 

a) dE/dt = Q.lAQrl{am)-'^rZl<i){i) 

b) dE/dt = l.lQrla-^/\mf/^mf'^\zf'^\{i) (13) 

c) dE/dt = l.&Qrla-^/\m)y^^mf'^\zT<f>{i). 

2.2. The star-star interactions and the inflow of stars 

In spherically symmetrical conditions, the influence of the 
central massive black hole on the stellar cluster structure 
and the inflow of stars into the MBH was considered by 
Peebles (1972) and solved by Bahcall & Wolf (1976) (here- 
after BW). BW showed that the presence of a MBH in a 
dense stellar cluster typically leads to the creation of a cen- 
tral stellar cusp with a density profile n{R) — no{R/ R^Z)'^ , 
with q ~ -7/4 and i?a - RssMbh/Mss- The stellar 
distribution function depends on the stellar energy as 
f{E) - E^/^, where E = GAfsH/r - l/2V^ and Vs is the 
stellar velocity. Both the creation of the cusp and the in- 
flow of the stars into the MBH are due to star-star gravita- 
tional interaction. The result was obtained for an isotropic 
stellar distribution function, but BW showed that it is 
valid if we use anisotropic distribution functions as well. 
However, in the case of AGNs with accretion disks, the 
interactions of the stars with the disk in the central part 
of the stellar system are inevitable, and the influence of 
the interaction on the stellar dynamics have to be taken 
into account. 

Let us compare the rate of the energy changes of a 

typical star (let it be the solar-type star) for the cases of 
star-star (hereafter s-s) interaction and star-disk interac- 
tion. For the former case one has 

d£;("") /dt ~ AwG'^M^V-^ng In A, (14) 
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where V^r is the effective relative velocity of the stars, rig 
is the stellar density and In A is the Coulomb logarithm; 
for the latter case Eq. applies. 

It is essential to note that in the present work we will 
use some effective values of physical variables. In principle 
these variables show a certain distribution at each radius: 
stars at each radius have various inclinations, velocities 
relative to the disk etc. A complex numerical approach 
would be necessary in order to model the distributions of 
these quantities. Knowing the distributions we would be 
able to compute accurately the mean or median values of 
the parameters. However, such a program is far beyond 
the scope of the present paper. Instead, we reduce the 
dimension of the problem by introducing a priori some ef- 
fective values at each radius, with values following from 
processes we consider. It permits us to determine physi- 
cally reasonable radial trends; but it is obvious that the 
obtained results are true only in a qualitative sense. 

If the mass of the CSC is larger than the mass of MBH, 
the effective velocity of the stars in the outer parts of the 
CSC is close to the virial velocity Vs « {GMsc/ RscY^^ , 
where M^c and Rsc are the mass and radius of the CSC. 
In the inner part of the AD the velocity at the almost 
circular orbits is close to V; « (GM/RY^^, where M{R) = 
Mbh + Msc{R) + Mab{R) is the total mass within the 
radius R. 

Ranch (1995) has shown that close to a black 
hole the characteristic two-body relaxation time tgs ^ 
v'^/{7rG^M^ns) 10-^^(ns/10'' pc~^) is much longer than 
the characteristic time of the stellar orbit alignment with 
the disk, ti ~ torhMg/ {ttR'^) /^d- This is still true in our 
case of more dense CSC if it is spherically symmetric, 
which can be seen from the numerical estimates of the 
correspondent energy-loss rates (in erg s^^): 

dE^'''>/dt ~ 5.9- 1032m^m8^^Vi4^nio(lnA)i, 
dE'^'^'^^ldt ~ 5 • 10=^5gS4?•|m^/^r:^/^ (15) 

where mg = M{R)/l(f Mq, nio = ns/(10^° pc'^). 

One can see that the star-disk interaction is relatively 
strong even in the outer parts of the AD. We can estimate 
the critical outer radius of the region where the two inter- 
actions became comparable, if the inclination of the orbits 
is negligible [sin(ief) ~ 1] in this outer region: 

r^J^' ~ 9.5i:'J\^J^m~'ml^^[nwiinA)i]-'/\ (16) 

that is i?"" ~ 0.1 pc (if j:^ ^ nis ^ mo = 1). 

At smaller radii the inelastic star-disk collisions be- 
come important, and the stellar orbits incline to the disk. 
If the orbit evolution starts from a non-inclined orbit, the 
orbit will quickly decline to the disk plane (as was shown 
by Ranch 1995). But when the inclination angle becomes 
small, the situation can change due to the increased stel- 
lar density. The reason is that the inclination of the orbits 
influences the relative velocities, and consequently the re- 
laxation times of the star-disk and star-star collisions. 



We will assume that the inclination angles of the orbits 
are smoothly distributed above the disk, and the distri- 
bution can be characterized by some effective inclination 
angle id- The characteristic relative velocities of stars in 
the co-rotating (with the disk) coordinate system are then 
close to VIsr — siniof, that is close to the average veloc- 
ity component normal to the disk. Introducing the effec- 
tive inclination angle into the expressions of the energy 
changes due to the star-star and star-disk interactions, we 
have from Eqs. (|^) and (^^ numerically (in erg s~^): 

dS'-^'^^/dt ~ 5.9- lO^^m^mg ^/^ri'2^nio(lnA)i/sin(icf),(17) 

dE^'-'^^/dt ~ 5 • lQ^'^QY.irlml'^rZl'^(t>{iod- (18) 

The important point is that the s-s interaction be- 
haves as l/sin(icf) and the s-d interaction is proportional 
to [sin(iof)]^, so both interactions may be comparable at 
small angles. It is intuitively clear, that an equilibrium 
situation exists when dE'^^^^ jdt ~ dE^^^^ jdt. The process 
is similar to the Brownian motion of the normal velocity 
component sini, and the equilibrium velocity dispersion 
corresponds to an equilibrium of braking and fluctuation 
forces in the Langevin equation. 

As the characteristic time of approaching the equilib- 
rium is relatively short (close to the characteristic time 
of orbit inclination, t\ ~ tdyn(-Ms/7r_Rg)/I]d), we will sup- 
pose further that the equilibrium condition dE'^^^^/dt ~ 
dE^'"^'>/dt is realized at any distance R < R"'\ This 
equilibrium means that in a dense enough CSC the stars 
will never be completely incorporated into the disk plane: 
when orbits are inclined too much, the pair stellar interac- 
tions will come into play, and they will scatter the stars to 
less inclined orbits (the eccentricities of the orbits will also 
be scattered to the same order of magnitude Cof ~ sini). 
The immediate consequence of this property, which leads 
to the fluctuation of the inclination angle around some ef- 
fective value, is that the fluctuating drag force maintains 
some constant dissipative flow of the stars to the center of 
the disk. The rate of the inflow is 

Vr = dR/dt = dR/dE{dE/dt), (19) 

where dE/dR = GMbrMs/R"^ erg cm" \ that is 
dR/dE ~ 3.37- 10^35^-1^-1^2 _^ erg^\ and dE/dt is 

the averaged dissipative energy loss of the star. We there- 
fore obtain 

Vr ~ 17Qi:4r^m~^ml^^rZl^^(l){i) cm s^^ (20) 

In stationary conditions, the flow of the stars is 

dus/dt ^ n ~ ATTnsVRR^ sm{ici), (21) 

where Ug is the stellar density at the AD plane and ic[{R) 
is the effective inclination angle of stellar orbits, which 
means also the part of the solid angle where the inflow is 
concentrated. 
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Numerically, this formula gives the flow: 



% ~ 2 • 10 ■^nioQS4r'g TTig ^m^^^r^^2 ^micf(f>{icf) yr ^(22) 

in units of stars/year, where nio is the stellar density in 
units 10^° pc"^. 



3. The distribution of stars in the central part of 
CSC 

3.1. Constant stellar masses and inflow of the stars 

As was shown in the previous section, the stellar orbits 
tend to decline to the disk and to shrink due to the energy 
loss of the stars crossing the disk, but the inclination of the 
orbits is limited by the scattering due to the gravitational 
star-star interaction. In the stationary case both processes 
are in equilibrium, which defines the inflow of the stars to 
the center. 

This process can be described with a system of two 
equations, which permit us to find the distribution of the 
stars. The first equation is the star inflow Eq. (p^, and 
the second is the condition of equilibrium of the inelastic 
(star-disk) and elastic (star-star) collisions. 



dE'^'''^^ /dt ^ dE'^''''^ /dt. 

Eq. ( p^ ) can be written in the form 



(23) 



sin(iof)(/i(ict)nio — 5 • 10''ny/{QT,4)msr^ ^m^ ^^'^r_l^^{2A) 

where hy is the stellar inflow in the units 1/year, uiq is 
the stellar density in 10^" pc^^. On the other hand, from 
Eq. (H), using Eqs. (|l^) and ([l|) one obtains: 



sin(ief)(/)(ief)/nio ~ 
1.2-10-3(QE4)-^ 



-2^3 



(InA)i 



(25) 



Eqs. ( p^ and (|2^) form a system of equations for two 
variables: nio(r) and ioi{r). Dividing Eq. (24) by ( p5| ) one 
has 



nio{r) ~ 0.6 • 103[rV/(ln A)l]l/2m^'m^/V:^/^ 



(26) 



This distribution is independent of the disk parame- 
ters, paradoxically at first glance. The reason is the un- 
derlying assumption of the equilibrium of non-elastic (s-d) 
and elastic (s-s) processes, but this condition is satisfied 
only in dense enough disks, so it is not valid for low surface 
density disks. 

Multiplying Eq. (p4|) by (psh, and using (for id < 7r/3 ) 



the approximation 0(i) = 2 sin icosi/{l -\- cosi) « sin i, 
one has 

sin(icf) « 

0.9[ny(lnA)i]i/6(QS4)-i/3(TOs/rs)'/'™8"'^''ri4" (27) 

This solution is valid in the innermost parts of the 
CSC/ AD system, where sin(icf) < 1; in the outermost 
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Fig. 1. Radial dependence of stellar density nio (continu- 
ous line) and sin(icf in the case of constant stellar flow for 
SS disk, constant surface density disk and NT disk (long- 
dashed, short dashed and dotted line, correspondingly 



parts the orbits are non-circular and we pose sin(ief) = 1 
there. 

If we use the disk surface density E of the NT disk 
from Rauch (1995), then 



sin(ief) ~ 0.14[ny(lnA)i]i/6(Q)-i/3x 
irn,/r,)y'm^'/'a%''{rn^r)-y''r'Jl 



(28) 



The dependences of sinicf on r for three zones of the 
SS accretion disk are: 

sin(jof)(r) ~ 

a)0.64Q-^^^{m,/rs)'^^Ha my^^ml^^^[hy{liiA)i]^/^rZl^^ 

&)0.34Q-l/3(TO3/r,)2/3(a)4/15(^)-l/5^-41/60^ 

K(lnA),]VV^/- 

c)0.3Q-l/3(^^/^^)2/3(^)4/15(^)-7/30^-ll/15^ 



[%(lnA); 



11/6^1/2 



(29) 



The dependence of sin(jof) and nio on r_2 for different 
E4(r) (the SS, NT and constant density disk cases) are 
shown in Fig. 1. We adopt Mbk = lO^M©, a = 0.01 and 
TO = 1. One can see that the stellar density rises quickly 
at small radii, while the effective inclination angle of the 
stellar orbits diminishes to the center. As can be seen from 
the picture, the difference in the behavior of sin(ief) for 
different disk models is not too large (though it has a wide 
minimum in the case of an SS disk), so further on we will 
use the NT model for our model calculations. 
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The inflow rate of the stars fiy in this consideration (in 
the approximate model we investigate here) is assumed to 
be given by a constant, which is determined by the con- 
ditions at the outer border, R = Rout- This outer border 
is close to the above-estimated value R'^"^ , determined by 
condition (|l^. Physically, the value n(i?out) is determined 
by the diffusion process, which depends on the proper- 
ties of the compact stellar cluster, and which for the case 
of a non-rotating CSC was described by Bahcall & Wolf 
(1976). It is clear that in the central part of the CSC the 
stellar density can be so large that physical collisions of 
the stars can occur; besides, the stars can lose mass when 
they cross the disk. These processes have to be taken into 
account. 

3.2. The contact stellar collisions 

Let us consider the change of stellar flow due to the 
contact (sometimes it is called 'physical') collisions of 
the stars; this is the non-elastic s-s process. The colli- 
sion velocities are about Vc ^ V sin i and the escape 
velocity from the surface of a star can be estimated as 

~ 6 • 10''[ms/rs]^/^ cm s~^. If Vc > w*, the averaged 
mass loss per collision is large (Murphy et al. 1991), and 
stars are completely disrupted if the impact parameter is 
less than some value. In the opposite case Vc < (if, say, 
sin(iof) << 1), the collisions can produce binary stars, or 
stars can be merged. Within the framework of the present 
paper we adopt equal masses of stars (which itself is far 
enough from reality!). In contact collisions a fraction of 
mass is always lost, but merging can increase stellar mass. 
Therefore, we will assume that both processes somehow 
compensate, and consider unchanged stellar masses, but 
take into account a decrease of the stellar density and 
the stellar inflow result of the direct s-s collisions. 

In Sect. 3.1 we have obtained formulae for the stellar 
density Ug and sin(icf), depending on the inflow of stars 
fly. These formulae are valid also for the case riy = hy{R). 
The frequency of the physical collisions of a given star 
with another stars (all stars have stellar radii Rs) is 

i/gg = AttUsR^Vs sin(ict), (30) 

and numerically: 

— 1 2 5/4 -11/4 . /. N _i /oi \ 

TOg r^m^ r_2 sin(Zof) s . (31) 

The frequency of the star-star physical (contact) collisions 
in the spherical layer is 

di^ss = I&t:^ RlnlV^^[fAii{ici)YdR, (32) 

thus giving, 

dv^s ~ 0.95 • lO^-^x 

rgm~^mgny(ln A)]^"'^rZ2 sin^(iof)c?r_2 ■ (33) 



The rate of the diminishing of the stellar inflow due to 
disruption of stars is d{ny)/dr = — gvdj^ss(l/s)(3-10^s/yr), 
where is equal to 2 or 1, for the disruptive or merg- 
ing case, depending on the Vr/w* relation. Besides, the 
cross-section of the stellar collisions depends on the rela- 
tive stellar velocity due to gravitational focusing at 
small velocities. We use these dependences in the numer- 
ical calculations (see below). Here we will neglect them 
and obtain approximate solutions with gv = 1- 

In the spherical layer between R and R + dR 

d{hy) ~ 

2.9 • 103%(ln A)^^TOs-2^2^2 sm'^{icf)rZldr^2- (34) 

Using sin(ict) from Eq. ([2^), one has {hy in units yr^^) 

d{ny) ~ 1.25 • lO^Ti^/^x 

{Q^i)-'/'{lnA)f/'m,-^/'r^^'ml^'rZl^^dr^2. (35) 

The solution of the differential equation can be ob- 
tained for different I](r) — Ar"^ . If we assume S = const, 
the solution is 

ny{R) ~ {rVy + 0.26 • W^{QJ:4y^/^x 

(In A)-/^m.-/Vf/3™^/«[.-/^ ^ (36) 
If we use S from the NT model, then 
d[ny{R)] ~ 0.5 • 102[7iy(i?)]^/3(Q)-2/3(ln A)7'/^x 

m72/3r2/3^5/6^8/^15(^_^^_7/15^-2^^_^^ (37) 

and the approximate solution is 

ny{R) ^ {h-;/^ + 3.4(Q)-2/3(ln A)-'/^ 

m.-2/V|/3mf al4'^(m)-Vi5[.:i - r,\]r\ (38) 

With these results and using Eqs. ( p6| ) and (p7|), we 
have the stellar density and the effective inclination angle 
of orbits 

nw{R) ~ 0.67- 103[ny(i?)/(lnA)i]i/2m^/V^/4TO7i, 
sin(icf)(^) ~ 0.93[ny(i?)(lnA)i]i/6x 

{Qj:,)-'/^m,/r,)y^m,'/'M/l (39) 

A more precise result can be obtained numerically if 
we take into account the dependence of the cross-section 
of the stellar collisions on the relative stellar velocity due 
to gravitational focusing at small velocities Vrei < Kssc — 
6 • 10''(ms/rs)"^/^ cm s^^. These cases take place at small 
sinief and/or at large distances, where we have to use 
collision cross sections as Css = 47ri?g(l-|- T4sc/K-ci)^- 

3.3. Stellar mass loss due to s-d interactions 

In the above considerations we neglected the change of 
the masses of stars due to s-s collisions, assuming that 
mass loss and stellar merging compensate and that stellar 
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Fig. 2. The radial dependence of the stellar density nio, 
and stellar flow riy in the case of Mbh = 10^ M© and 
the values of the initial stellar flow equal 0.1, 1.0 and 10 
stars/yr (dashed lines). 

masses keep their initial values. This approximation is rea- 
sonable at large radii since the stellar density is small, and 
collisions of stars are rare. On the other hand, stars would 
systematically lose mass when crossing the disk, and this 
process is more intense at small radii. 

The frequency of the star-disk interactions of a single 
star in a circular orbit with radius R is inversely propor- 
tional to the half-period, v^d = ^/^T{R); numerically 

^.i^^O-Tl-lO-V^/^rZ^/^s-i. (40) 

The differential frequency of the star-disk crosses 
by stars from the volume of a spherical layer dV = 
AttR"^ sm{ic{)dR is 

di^s^ = AnnsAT^^R^ sm{ici)dR, (41) 
and expressed in [s^^] gives 

di^sd ^ 0.6m~^ml^^[hy/(\iiA)i]^/^rZl^^sm{i^f)dr.2- (42) 

As the effective inclination angle is non-zero (and typically 
it is larger then several degrees), the star-disk velocities 
are supersonic and the stars lose their masses when cross- 
ing the disk. Let us suppose that the mass-loss rate is 
proportional to the dissipated energy 

dMs/dt = q,^{dEl/dt){dM/dE), (43) 

where dEl/dt is the rate of the energy dissipation in the 
disk and dM^/dE = -Rs/{GMs). 

Obviously, the envelopes of the red giants will be blown 
off in the process of crossing the outermost parts of a 



typical accretion disk of the bright AGN, so further on 
we will consider only the main-sequence stars. Denoting 
^ 10^^(7m-2, and using Eqs. (^ and (||) we have (in 
units of g s^-^) 

dM,/dt = -2.6-10i«ft„_2r-3^m-iQS4m^/'r:^/'</)i(ief).(44) 

We can deduce the radial dependence of nig from 
dMs/dR — {dMs/dt)/Vr. Assuming cos«of — 1) and 
Tg ~ iris (appropriate for the main-sequence stars), one 
has 

dms/ms ~ ~2.3qm-2'ms{dr^2/r'i2)- (45) 
The solution of the equation is 

rus ~ nisQ exp{-ii:[l/r_2 - ^/r'^-l]}- (46) 

where K — 2.3g,n-2^7i8 5 and r^^l = r™* is the outer (crit- 
ical) radius of the stellar inflow. The peculiarity of the 
process is that (within our approximations!) it does not 
depend on the S4(r) behavior. 

In the region r_2 << ^''^2 solution is rris ~ 
rriso exp{— _fir/r_2}, and one can see that the mass of a 
star quickly goes to zero when r_2 < K. Physically this 
means that at this distance the time scale of the mass loss 
becomes less than that of the radial drift (inflow) , so the 
stars with diminished masses crowd close to r_2 ~ K. The 
distance i?ni ~ O.OlKqtn-2 pc is the radius of melting of 
the stars due to the star-disk crossings. 

It can be noted that the time between consecutive 
crossings is several years at this distance (mg = 1), and 
the duration of crossing is about one day, so the stars re- 
ally can be bloated due to deformation and heating of the 
surface, and this leads to even quicker melting of the stars. 

Adopting a minimal stellar mass of ~ O.OIM©, one 
has (r_2)min ~ -K/\n{mso/'ms), that is (r_2)min ~ 
0.5g„i_2»7i8 for m^o = 1, and (r_2)min ^ qm-2ms for 
TOsO = 0.1 . 

As the stellar density drops quickly due to the melting 
of stars at r_2 < (?'-2)min, in the mmierical calculation 
we used the approximation 

nio ~ 1.3- 103[riy/(lnA)i]i/2x 

ml^\zT^7o exp{[if/r_2(l - 0.5K/r^2]}, (47) 

that is the stellar density is enhanced close to r_2 ~ 
2.3gni-2?7i8 and then drops at smaller radii. This means 
the appearance of some 'stellar hoop' with the radius 
r-2 — 2. 3(7111 -2 "^8- The numerical calculations of the emis- 
sion line profiles (see Sect. ^) give reasonable results with 
(7ni-2 — 1, which we accept below as the standard value. 

The influence of the process of star melting at small 
radii was taken into account in the numerical calculations 
of the models, and the results are presented in Figs. 2 and 
3. 

Our approximation works for sufficiently dense disks, 
so we show the results of the numerical calculations for 
the case of NT disk with a = 0.01 and m = 1. 
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Fig. 3. The radial dependence of the stellar density nio, sin(iof) and stellar flow riy for different black hole masses: 
a)MBH = 10'' Mq, b)MBH = IO^Mq, and c)Mbh = IO^Mq. Thin continuous line shows the nio distribution in the 
absence of stellar collision effect. 



In Fig. 2 we show the radial dependences of stellar 
densities and stellar flows for the case of Mbh = lO^M© 
with different initial (at the outer border) values of the 
stellar flow hyQ = 0.1, 1, 10 yr~^. One can see that at small 
distances all the flows ny{R) become almost the same and 
independent on the initial values of riyo at large radius. 
This result is notable, as it shows that in fact the accretion 
disk itself can be built with the gas resulting from the s-s 
disruptive collisions in the dense enough CSC. 

We used this hypothesis for calculation of the refer- 
ence models of different AGNs, assuming that the accre- 
tion rate onto the MBH is always equal to the initial mass 
inflow by stars at the outer border of the CSC/ AD system, 
hy-ms = m^rh. In Fig. 3 we show the radial dependences of 
stellar densities, inclination angles and stellar flows for the 
reference models with Mbh = IO'^Mq, 10®Mq, IO^Mq; 
also cases without contact stellar collisions and mass loss 
due to the disk crossing are shown for comparison (thin 
lines), and the behavior of the stellar inflow fiy{R) is 
marked (dashed lines). 

4. Broad emission line region 

There is still no commonly accepted physical idea about 
the nature of the broad emission line regions (BELR) of 
AGNs. The possible role of the accretion disk was stressed 
by Dumont & CoUin-Souffrin (1990) and it was subse- 
quently considered in a number of papers (e.g. Collin & 
Hure 2001). There were, however, also many attempts to 
involve stars (viz. excited, bloated or comet-like stars) in 
the explanation of problems of clouds, moving and emit- 
ting in BELRs (Edvards 1980; ShuU 1983; Penston 1988; 
Scoville & Norman 1988, 1995; Kazans 1989; Zmek et al. 
1994). The most recent was the study by Alexander & 
Netzer (1994, 1997), where the authors considered emis- 
sion from the outer envelopes of excited stars and used a 
spherically symmetric model of the structure and dynam- 
ics of the CSC in the vicinity of the MBH (Murphy et al. 



1991). The application of a spherical model to the stellar 
dynamics in BELRs is questionable, because it completely 
ignores the interaction of stars with the AD. Our model of 
stellar dynamics includes both interactions which together 
determine the correct spatial and kinematical distribution 
of stars and we can now built a model of the BELR based 
on that knowledge. 

To define the physical nature of the BELR clouds, and 
the model of the star-disk interaction in more details, we 
refer to the paper by Zurek et al. (1994), where the struc- 
ture of the plasma around the star moving through the 
disk is presented as a combination of a hot bow-shock 
filled with hot plasma, and a more dense and cold wake, 
or tail, behind the star. We suppose that when the star 
comes out from the disk surface, the appearance of the 
bow-shock from the surface of the disk can be treated as 
a fireball, and the stellar wakes can be identified with the 
clouds of the BELR. 

Zurek et al. (1994) assumed that the BELR clouds are 
actually wakes of matter, dragged by the stars from the 
disk. We accept this idea in a slightly changed variant, 
taking the view that the matter dragged from the disk 
creates a fireball, and the wakes contain the mass, lost by 
the stars when they cross the disk (both parts are the same 
to within an order of magnitude) . As was described in the 
previous section, the fraction of mass lost by the star is 
proportional to the dissipated energy AE, that is AMg = 
0,01qm-2AEdM/dE, where dMjdE = R^/{GM^). 

Taking AE from Eq. (|), one has 

Aft ~ AM, ~ 3.5 • lO26g^_2QS4r3m-im8rZi0i g, (48) 

where q-m-2 = 1 means that we assume that 1% of the 
star's energy dissipated in the disk is used for the mass 
loss of the star. 

As the balance of forces acting on the wakes depends 
on the parameters of AGNs and is model-dependent, we 
assume here the simplest case: the wakes follow the stars. 
Let the wakes be cylinders with lengths approximately 
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equal to the height of the stellar orbit above the disk, 
It ^ i?sin(zof), and let the radii of the cylinders be close 
to the stellar radius, rt ^ TsRq (this last suggestion may 
require the presence of a hot corona to provide the cold 
gas confinement). Then the mean particle density in the 
cylinder is 

n ~ Mt/[ltTTr'^ sin(ief)mp], (49) 
or numerically 

n ~ 3.8 • 10^^(7_2QS4rsTOjr^?7i8rZ2</'i/(sinicf) cm^"^, (50) 
and the column density is 

TV ~ 6.5 • lO22g_2<3S4r,m7im8rZ^0i/(sinw) cm~^ (51) 

The covering factor (CF) of a single cloud (along the 
disk radius) is 

n = 2i?Jt/47ri?2 ~ 3.7 • 10~VsSiniefA-2, (52) 

and the total covering factor of many clouds in a spherical 
layer from R to R + dR (the differential CF) is 

CF = AnnrisR^ sin i^fdR, (53) 

where Ug is the stellar density (the number of clouds per 
unit volume is close to the stellar number density because 
we assume the life-time of a cloud is close to half the or- 
bital period of the star: one star - one cloud). Taking 
into account that at small and intermediate inclination 
angles 0i/(siniof) ^ siniof, and that the flow of the ioniz- 
ing photons is proportional to <I>i ^ 10^'^msmrZ2y one can 
see that the particle densities and the radiation flows are 
within the range of calculations of Korista et al. (1997) 
for the atlas of equivalent widths {EW) of quasar broad 
emission lines. These calculations give EW as a function 
of two parameters, $i and n, supposing that all the clouds 
have column densities N — 10^^ cm~^ and the total cov- 
ering factor of the clouds is CF = 1. This permits us to 
calculate the profile of an emission line as 

I{V) K EW{^u n)VlngR^ sini^fdR 

d{V - Vp- Vch)}d(pdvch (54) 

where V is the velocity projection along the line of sight, 
X is a constant for a specific fine, EW{^i,n) is the EW 
of the line calculated in Korista et al. (1997) as a function 
of $i and n, is the coverage factor of a cloud at the 
distance R, Ug is the number density of stars (and so of 
the clouds), Cr is the relativistic correction, Vp{ip,9) = 

cos (/? sin is the projection of the rotational velocity 
onto the line of sight, {ip is the azimuth angle and 6 is the 
angle between the line of sight and the axis of the disk), 
Vch is the projection component of the chaotic velocity of 



stars with dispersion D{R) ~ 0.67 • 10^(Af8/r2) sin^cf cm 
s~^, where M{R) is the total mass within radius R and 
(5() is the delta-function. 

Note that the velocity dispersion D of the chaotic mo- 
tion contains sinief, because the velocity dispersion is due 
to the distribution of the inclination angles and eccentric- 
ities around the effective values; we assume here that the 
distribution of the projection velocities is Gaussian in the 
coordinate system co-rotating with the AD. 

We convolve the rotational velocity with the chaotic 
velocities of the stars; this means that the velocity of 
clouds contains both the rotational (disk) and chaotic 
velocity components of the stars. i?min and i?max in the 
integral hmits are taken correspondingly as O.OOlpc and 
(mgm) ^/^pc, the last is approximately determined by the 
inner radius of the putative obscuring torus. 

We performed the computations for the disk surface 
density profile of NT given by Eq. (0). 

The results of the numerical calculations of the pro- 
files of the Ha line, normalized to 1 at the maximum, are 
presented in Fig. 4. Each panel contains three profiles, 
corresponding to different inclination angles of the line of 
sight to the axis of the AD: 9 = 10°, 30° and 60°. The 
general shapes of the profiles are determined by the com- 
bination of a broad part, which depends mostly on the 
rotation of the inner region of the disk, and a more nar- 
row central part, which is determined mostly by stellar 
chaotic motions in the outer parts of the system. 

In Fig. 4(a-c) we show line profiles for the Eddington 
mass inflow {m — 1) and for different masses of the MBH, 
rrig = 0.1, 1, 10, which would correspond to a Narrow-Line 
Seyfert 1 galaxy (NLSl), a Seyfert galaxy and a quasar, 
respectively. Our calculations support the idea that the 
narrow line Seyfert 1 galaxies are those with small MBH 
mass and almost Eddington mass inflow; or they can be 
objects with intermediate MBH mass and with a mass 
inflow less than the Eddington one, seen at an angle 9 < 
30°. 

Our theoretical consideration of the stellar distribution 
above the AD gives reasonable line profiles, supporting the 
idea that the emitting clouds are somehow related to the 
stars. They can be the wakes, following the stars after 
crossing the disk, or/and some gas envelopes of the stars 
themselves, which are disturbed by disk-crossing and by 
close passages of other stars. The mean distance between 
the stars at stellar density rig ^ 10^ pc~'^ is dr ~ 10"'^ pc, 
and many stars pass at even smaller distances, exciting 
each other by tidal forces and tangential collisions. How- 
ever, in the present paper we used the specific model of 
clouds as the stellar wakes. 

In our model the covering factor of the clouds along the 
AD plane became higher than 1 at large radii, so the spec- 
trum of ionizing radiation has to be changed in the outer 
parts of the BELR. Of course, our model calculations are 
preliminary and rough, and a more detailed model has to 
be considered for the sake of comparisons with the ob- 
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Fig. 4. Ha line profiles, calculated with Mbh = lO'^MQ(a), Mbh = lO^M0(b), and Mbh = IO^Mq(c). Dotted, solid 
and dashed lines correspond to ^ = 10°, 30° and 60° 



served profiles. But the common property following from 
the physical sense of our approach, viz. that the wings 
of the profiles should be more variable than the central 
peaks, seems to be confirmed by reverberation investiga- 
tions. 



5. Variability 

In the present section we will estimate the level of variabil- 
ity driven by the stars. Obviously, the interpretation of the 
variability demands the calculation of rather complicated 
models, but here we will make a simple analysis of the 
variability expected from the inelastic processes of star- 
disk and star-star interactions. We consider two types of 
processes, connected with the energetic events of s-d and 
s-s collisions: flashes of UV radiation and X-ray flares. 
We consider a very simple model: let every event produce 
a fireball, filled with matter dragged from the disk (s-d 
crossing) or from stars (s-s collision), and expanding in a 
hot coronal medium. When the fireball becomes Compton- 
transparent, it produces a UV flare. The radiation of the 
flare interacts with the hot electrons of the corona above 
the disk, thus producing a flare of the up-scattered X- 
radiation. 

The real picture is of course much more complicated. 

The fireball radiation spectrum evolves; quickly expand- 
ing fireballs create blast waves in the environment, and 
the bow shock can be accelerated when propagating from 
the disk surface into the corona with diminishing den- 
sity. It is possible also that the matter of the fireball is 
heated later on to the coronal temperature, thus creating 
a transient patchy corona. The undisturbed, steady-state 
corona (without any influence of stars) can extend up to 
R = 300i?g (Kawaguchi et al. 2001); the patchy corona 
created by s-d interactions can extend to even larger dis- 
tances. 

All these processes require a more detailed description, 
but here we consider the simplest one, based on the simple 
picture of fireballs. 



5.1. The star- disk events 

When a star comes out of the disk, it is preceded by a 
bow-shock, accompanied by a fountain of the matter from 
the disk. The total amount of gas thrown out above the 
disk consist of two parts: AMd is the mass dragged out 
from disk by the star, and AMg is the mass lost by the 
star. 

We assume that the fireball is created by the first mass, 
and that the second mass creates the wake, following the 
star and showing up as a BELR cloud, as was discussed 
in the prc;vious section. 

Then the mass of the fireball is Mf = AMd = 
gi7ri?gSd, where < 1; numerically 



Mf ~ 1.6 • IQ^^qirlY^i g. 



(55) 



Let the fireball contain the energy deposited by the 
star into the disk matter, 



A£; ~ 0.7 • lQ^^Qrli:AmiirZl(t>{i)/{cofii) erg 



(56) 



We assume that the fireball expands and produces the 
flare when it becomes transparent, that is optically thin 
relative to the electron scattering. At this moment its ra- 
dius fits the relation Ripa = 1, where p = M{/{4:/3ttR^) 
is the density of the fireball, and a ~ 0.4 cm^ g~^ is the 
Thomson scattering cross section per gram of the ionized 
cosmic plasma. Then the radius of the fireball at the mo- 
ment of the flare is 



Rf ~ (MfCrT/4)^/^ ~ .4 • 1013(^^5^^)1/2^^ 

The density at the moment of the flare is 
Pf ~ Mf/(4/37riif 
so the particle density is 



6.3•10-lV'^V'S-^/^ 



3.2- 10"gi 



l/2_^_l„-l/2 



(57) 



(58) 



(59) 



and the column density is close to Nn = 10^^ cm 
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The temperature in the fireball changes due to its ex- 
pansion approximately as 



Tf ~ [giAi;/(4/37rai?''')]i/4 K 



(60) 



where a = 7.6 • lO^^'^ erg cm^"^ K^^, and R is the radius 
of the expanding fireball. 

At the moment of the flare, the temperature is there- 
fore 

Tf ~ .8 • lO^Qxr-iS- ^/'m^/V^/'0li/4 K, (61) 

where Qt = qi~^^^Q'^^^, and the spectral maximum of 
the flare is in the far UV-band, 



1/ - 5 • 10'^ Hz. 



(62) 



The estimate of the minimum duration of the flare 
(assuming that the radiation leaves the fireball when the 
latter becomes Compton-thin) is 



Atf < 2Ri/c ~ .3 • lO^^iS^ s 



(63) 



and the average luminosity during the flare is 

Lf = AEf/Atf < 3. • lO^igr.Sy^TOsrZ^?!)! erg (64) 

where q = qiQ. 

In the NT disk we have 

Lf ~ 0.5-10*^qr,ml/\zT{m,)y^°rZl'^Ui erg s-i.(65) 

The real duration of the flare is rather longer, and 
the luminosity is less than in these estimates. When the 
temperature of the fireball becomes lower than 10^ K, the 
plasma in the fireball partly recombines, which increases 
the opacity to above its Thomson value. Besides, if the 
fireball expands with a velocity close to the sound velocity 
Vf ~ (fcB/m-pTf)^/^ ~ 3-10^ cm s~^, the time of expansion 
is tcx — Ri/Vf — 10® s. So, rough estimates show that the 
characteristic duration of the flare is in between the two 
limits 



3 • lO^^iEyVs < Aif (s) < 10®qiSy^-,/(Tf5) 



\l/2 



(66) 



Now we estimate the combined effect of all star-disk 
events. 

It has to be noted here that the results of model calcu- 
lations of luminosity and variability, connected with the 
star-disk interaction, depend strongly on the values of the 
rates of energy change, estimated by Eqs. (|lj) and (^. In 
principle, these estimates can be increased, because the co- 
efficient Q in dE'-^'^^dt is unspecified, and dE^^^^/dt can 
be increased, e.g. due to collective processes in the flat 
distribution of stars. The same effect of increasing lumi- 
nosity is expected if the equilibrium dE'-^'^^ /dt ~ dE^^^^ /dt 
is not reached, due to larger values of sin(iof). Because of 
this, we calculated all quantities of interest (see below) for 
our reference models with wg = 0.1, 1, 10 for two cases: 1) 
energy change rates given by Eqs. (|lj) and (||); 2) energy 
change rates increased by a factor of 5. 



The mean luminosity provided by star-disk interaction 
in the ring between R and R + dR is dL = AEfdvgd, with 
dvsd from Eq. (E2): 



2 - 10'*^QE4K/(lnA)i] 



1/2, 



r~mg"mg' r_2 ' </'(*) sin(icf)rf''-2, 
and the total luminosity of the flares is 



LfT = / LfAifdi^sd = / AEfdv, 



(67) 



(68) 



From calculations of the reference models with 7713 = 
0.1, 1, 10, we have the total UV-flare luminosities (in units 
of 10^3 ergs"^) in the first case Ln ^ 0.02,0.07,0.39, 
and for the second case {dE/dt 5 times larger): Lfr — 
2.2,29,7.4. 

One can see that in the second case the total luminosity 
of fiares resulting in the s-d interactions is comparable to 
the variable fraction of the UV luminosity. 

The average number of the simultaneous flares in every 
dr -ring is dL/Lf, so the total number of coexisting flares 
is 



N = Atidv, 



sd- 



(69) 



For the three MBH masses (wg = 0.1, 1, 10), we have (with 
Ai'"™) in the first case ~ 3.6,5.7, 100, and in the sec- 
ond case A^ ~ 17, 73, 300. 

The weighted luminosity of a flare is Ff — Lfr/N, and 
the mean-square deviation of the luminosity due to chaotic 
variability is 



DL 



sd 



(70) 



Corresponding values for the three MBH masses are in 
the first case DL/Li = XjN^l'^ ~ 0.5,0.4,0.1, and in the 
second case DL/Lf ~ 0.24,0.12,0.06; so, the variability 
amplitude diminishes with increasing luminosity of AGNs. 

5.2. The star- star events 

Similarly we can estimate the variability connected with 
the star-star contact collisions. We assume that the 
average mass-loss per collision is 
Mf ~ O.lC-iAfs = 2 • 1032C_ims g. 

The parameters of the s-s fireballs (the total energy, 
the radius of the fireball, the temperature of UV radiation, 
the minimum and maximum duration of the UV flare are: 

AE = 0.45 • 10^"C-i"ism8rZ^sin2(iet) erg, 

i?f ~ (MffTT/4)l/2 ^ i ^ . 10l6„^l/2^1/2 

Tf ~ 3.6 • lQ\^J^m~^/^ml'\-^/^ sin^ i,, K, (71) 

Atf" > 2Rf/c = 0.7 • iovy2(^yi^ s, 

Atr" < 2i?f/14ound ^ 3 • 10^ s. 
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The corresponding average luminosity is given by 
LY'^^ < .64 • IQ'^^Q'J^ml/^msr^^ sin^ ief erg s-\ (72) 

Lf " > 2 • lQ^\]!^mll'^ms,r-^ sin^ iefergs-^ (73) 

By definition, the column density of the fireball at the 
moment of the flare is = Rn ^ 10^"' cm~^, so n ~ 
10^ cm^^. 

Taking dv^s from Eq. (|3^) , we have a number of simul- 
taneous flares N — J Aifdiyss, which is typically less than 
1. This means that s-s events lead to even lower average 
luminosity than s-d ones, though the amplitude of the s-s 
flares are larger. 

5.3. The relation of the UV and X-ray events and the 
frequency dependence 

Of course, the variability connected with stars may not be 
the only mechanism operating in AGN. Their variability 
can be mostly driven by the variable (unstable) accretion 
rate onto the MBH. The process of accretion itself is rather 
complicated, and the theory develops quickly (see the re- 
cent review in Park & Ostriker, 2001). The variability ob- 
servations show some common properties of the variabil- 
ity of the AGNs and the stellar-mass Galactic black hole 
candidates (Utley & McHardy 2001, Chiang et aU. 2000), 
and for that class of objects the model of distributed flares 
has been considered (Galeev et al. 1979, Czerny & Lehto 
1997, Poutanen & Fabian 1999). On the other hand, it 
has been shown that UV and X-ray variability of AGNs 
cannot be explained with models assuming simple repro- 
cessing of the X-ray emission of a variable point source 
placed at the symmetry axis above AD. Is it possible that 
s-d collisions can work in AGNs instead of the magnetic 
flares supposed to exist in the Galactic black holes? 

From the rough estimates provided in the previous sub- 
section one can conclude that the s-d interaction can lead 
to a noticeable variability amplitude only if both rates of 
energy dissipation - due to s-d and s-s interactions - are 
enhanced several times in comparison to our basic esti- 
mates. This can possibly be achieved with a more realis- 
tic distribution of stellar masses. An increase might also 
come due to unusual properties of the stars in the very 
compact stellar clusters. Deviations from the equilibrium 
conditions (which was assumed by us for the approximate 
solutions) also promote larger values of sin(ief) and, as a 
consequence, lead to more power radiated from s-d inter- 
actions. 

Obviously, the UV events are divided into two types: 
those connected with s-d and with s-s interactions. The re- 
lations of the two types of interactions to the X-ray events 
are different. 

In the s-d case, the UV flares occur at relatively small 
heights above the AD, which are close to the fireball size, 
h ~ 0.39 • 10^^(giS4)^/^rs cm. The time for reaching this 



distance is t ~ R{/V{, where R{ and V{ are the size of 
the fireball and the velocity of its expansion, and the time 
from the appearance of the blast wave to the maximum of 
the UV luminosity is t ~ 0.39 • 10'^{qiT,4)^/'^rs s. The size 
of the region of up-scattering is about the coronal scale 
height, around 10^^ — 10^^ cm, which leads to a time-lag 
between UV and X-ray flares of about Atf ^ 10^ — 10^ s. 

In the s-s case, most events occur at large heights above 
the disk, so the fraction of photons reflected by the corona 
is small, and the duration of Compton-scattered X-ray 
events, as well as the time-lag, are much longer than in 
the case of the s-d event. 

Using the estimates from the previous section, we 
can consider the amplitude-frequency relations following 
from star-disk and star-star dissipative interactions. In 
Fig. 5 we plot the dependences of the UV flare luminosity 
if = AE{{R)/Atf{R), divided by frequency dv{R) on the 
frequency. We plot both Lfss/di'ss and Lfsd/di'sd in the 
same figure using the above-defined dependence of both 
values on radius and replacing the radial coordinate with 
the frequency coordinate where the frequency-radius rela- 
tion is taken from the determined characteristic frequency 
of events at a given radius. Let us call the dependence 
"frequency-amplitude dependence" (FAD). Note that the 
FAD cannot be strictly compared to the power density 
spectrum (PDS), observed in the UV and X-ray bands 
(e.g. McHardy & Czerny 1987, Chiang et all. 2000). These 
PDS are obtained by Fourier transformation of the light 
curves, which are usually measured in frequency bands 
about 10""^ — 10~^ Hz and show a sort of red noise, that 
is PDS oc z^~P with p ~ 1 — 2; sometimes p ^ is found 
at the low-frequency limit. The FAD spectrum shown in 
Fig. 5 has a different meaning: it is just the dependence 
of the flare luminosity per frequency interval on the fre- 
quency of the events. Besides, we recall that we use the 
equal stellar mass and other approximations, so the real 
slopes of the curves could be steeper and the curves could 
cover a broader frequency range. Nevertheless, some gen- 
eral prediction about PDS can be derived from our calcu- 
lations: 

i) the PDS of AGNs can be divided into two parts: 
fast variability (on time scales from hundreds of seconds 
to about a year, which is connected with the s-d processes) 
and slow variability (on time scales from several years to 
thousands of years, connected with the s-s collisions). 

ii) the PDS slope can change abruptly somewhere in 
between the fast and slow part, and the slow variability 
part of the PDS has a larger amplitude. This conclusion 
seems to be in agreement with the results of long-time 
observations in X-ray band (Markowitz & Edelson 2001), 
and it is definitely supported by long-time observations in 
the optical band (Czerny et al. 1999, Czerny et al. 2001). 

iii) the differences of PDS between the slow and fast 
parts may depend on the AGN luminosity and other pa- 
rameters. 
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Fig. 5. The Frequency- Amplitude dependence, including 
s-d (higli-frequency) and s-s (low frequency) interactions 

Our preliminary conclusions are that the presence of 
a CSC in AGNs docs not contradict the observations, the 
CSC can contribute to AGN variability in some cases, and 
further comparison of variability of AGNs and Galactic 
black holes is needed. 

6. Discussion and conclusions 

The role of the compact stellar cluster in the creation and 
evolution of AGN is still poorly understood. We suppose 
that this role is very important, so that the AGN phe- 
nomenon can be considered as a result of interactions of 
three main physical subsystems: the CSC, MBH and gas 
subsystem. We call this approach the interacting subsys- 
tems theory (Vilkoviskij et al. 1999). 

Our model of the CSC-disk interaction, though very 
approximate, is the first one which treats star-star and 
star-disk interactions jointly, so it can serve as a starting 
point for considering some new physics in AGN problems. 

We showed that the star-disk interaction tends to order 
stellar orbits, dragging them to the disk plane and making 
them circular, but the star-star gravitational interaction 
scatters these orbits and prevents them from being incor- 
porated into the disk. We assume equilibrium of the two 
processes and use a simple inflow model instead of analyz- 
ing the kinetic equations. The other simplification was the 
assumption of the equal stellar masses, which is acceptable 
as a first step only. 

Our results can be compared with the results of Rauch 
(1995, 1999). In the last work he calculated numerically 
the evolution of the spherically-symmetric stellar cluster 
with stellar collisions, but without s-d interactions. He 
finds that the internal stellar density distribution achieves 



a quasi-stationary state in the late stages. Our calculations 
of the radial stellar density distribution with direct stellar 
collisions and s-d interactions are in qualitative agreement 
with those of Rauch (1999). 

On the other hand, our calculations of the shape of the 
stellar density distributions in the central parts of AGNs, 
are naturally different to the results of Rauch (1995), 
where he considered the evolution of individual stellar or- 
bits crossing the disk (without s-s interaction). He calcu- 
lated the evolution of the stellar system, which surrounds 
the AD in the first moment, and showed that the evo- 
lution leads to the completely flat system, incorporated 
into the disk plane in the final stage. In contrast, we take 
into account (besides the s-d interaction) both elastic and 
non-elastic star-star interactions. This leads to a quasi- 
stationary distribution of stars, which is flat only in the 
central part of the disk. 

The determination of a more realistic stellar distribu- 
tion in the innermost part of the central cluster allowed 
us to consider some consequences of the complex system 
of AGN, containing MBH, AD and CSC. 

We showed that the simple model of stellar wakes ac- 
companying the s-d interactions gives profiles of broad 
emission lines similar to those observed, providing a rea- 
sonable solution for the puzzles of both generation and 
dynamics of the BELR clouds. We used even further sim- 
plified models for estimates of the variability connected 
with s-d and s-s interactions, and we come to the pre- 
liminary conclusion that a noticeable role of the stars in 
the variability is possible, but only if both dissipative and 
scattering interactions of stars are enhanced relative to our 
preliminary estimates. The variability produced by stars 
is naturally divided into a fast and a slow parts, corre- 
sponding to s-d and s-s processes. 

Of course, the main questions remains whether the 
stellar density of CSC in AGNs is large enough for all the 
effects discussed to proceed efficiently. Here we assiimed 
that the masses of CSC in the bright AGNs are about an 
order of magnitude larger than the MBH mass. 

The real stellar density and masses of the CSC in 
AGNs are still unknown. In the paper of Schinnerer et al. 
(2001) the near-infrared high angular resolution speckle 
imaging of the central region of the Seyfert galaxy NGC 
3227 was performed. The authors showed that the CSC do 
exist in the inner 300 pc, and that in the innermost (unre- 
solved, < 70 pc) part there is the even more dense core of 
the stellar cluster, which contributes less to the luminos- 
ity, but contains 85% of the mass of the cluster. This result 
supports our leading idea of large enough masses of CSCs 
in AGNs. A fuller investigation of 112 Seyfert galaxies in 
1.6/im with the Hubble Space Telescope was performed 
by Quillen et al. (2001). They showed that there are un- 
resolved sources in all Seyfert galaxies with luminosities 
i(1.6 iim) ~ (10^1 - 10^2) erg s'^ correlating with 0[III] 
A5007 and with hard X-ray luminosities for Seyfert 1.0- 
1.9 sources. The luminosities are much lower in Seyfert 
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2.0 sources. Due to this trends the authors concluded that 
the unresolved sources are mostly of non-stellar (hot dust) 
origin; but we note that in our model the correlation of 
the AGN luminosity with the mass of the CSC is ex- 
pected. More detailed observations are needed for solving 
the question about CSC masses in AGNs. 

It follows from our models that the accretion disk it- 
self can be created as the result of mass loss due to con- 
tact stellar collisions in the rotating, compact and mas- 
sive CSC. In this case the collision disruptions of stars can 
partly feed the accretion disk flow, and another part of this 
gas (the hot part) can support the outflow of the hot gas 
from AGNs, as was supposed by Vilkoviskij et al. (1999). 
This opens new possibilities for the more self-consistent 
models of the AGN structure, where the accretion disk, 
obscuring torus and polar outflow are created (partly at 
least) due to stellar collisions in the rotating CSC. The 
transfer of angular momentum from AD to the CSC also 
has to be taken into account. 

One can imagine the following general picture of the 
AGN's duty cycle: 

The bright phase of AGNs should be relatively short, 
about several times 10^ years, and usually the activity 
cycle consists of three phases. It starts with creating a 
new CSC in the dense gas cocoon around the MBH (say, 
generated by a galactic merger), which corresponds to the 
IR-bright IRAS galaxies (the first phase) . In this phase a 
new CSC is formed due to the powerful star-burst in the 
galactic center, which leads to the creation of a new AD, 
thus initiating the star-disk dynamics, similar to what was 
considered by Ranch (1995). 

The subsequent growth of the AD luminosity, the cre- 
ation of the central flat stellar distribution, and the ap- 
pearance of the polar jets and the hot-gas outflows make 
polar holes in the cocoon, and the second phase of activ- 
ity sets in; the AGN looks like a typical QSO or Seyfert 
galaxy with an obscuring torus. Now in this phase most of 
the stars in the polar regions of the rotating CSC are ex- 
hausted, and stars cross mostly the periphery of the AD 
due to the diffusion into an effective loss-cone. The AD 
itself can be supported by stellar collisions in the rotat- 
ing CSC, and our solutions approximate stellar dynamics 
in this phase. During approximately 10^ years the mass 
of the CSC and the stellar inflow will gradually decrease, 
thus leading to the third phase with diminished activity, 
corresponding to the weak Seyferts and LINERs (galaxies 
with Low Ionization Narrow Emission Region) . 

Therefore, the results of our paper correspond to the 
second phase, the 'adult and strong' age of an AGN. 

It is clear that more detailed model calculations are 
required. Our description of star-star interaction in the 
non-spherical, flat stellar system is very rough and the 
results should be verified with n-body simulations. Ad- 
vanced models should also take into account the distribu- 
tion of stellar masses and a more accurate description of 
the accretion disk. New achievements in theory and obser- 



vations of AGN show that the accretion flow is possibly a 
combination of a geometrically thin cold disk in the outer 
part, and an advection-dominated accretion flow (ADAF) 
in the central part of the disk in some cases (Ichimaru 
1977; Rees et al. 1982; Narayan & Yee 1994; Abramowicz 
et al. 1995, Loska & Czerny 1997), with the hot corona 
still present above the cold disk (Esin et al. 1997; Esin et 
al. 1998, Witt et al. 1997). 

However, even the approximate results presented here 
show that compact and massive stellar clusters in the cen- 
ters of AGNs can play an important role in the creation 
and physical properties of the AGN engine. 
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